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TECHNICAL MEMORANDUM 1319 


TORSION AND BENDING OF PRISMATIC RODS OF 
HOLLOW RECTANGULAR SECTION* 

By B. L. Abramyan 

In the present paper a solution is given for the problem of the 
torsion and bending of prismatic rods of hollow rectangular sections. 

As in the former paper (reference l) , the method given by 
N. Kh. Arutyunyan (reference 2) of introducing auxiliary functions was 
employed in the solution of this problem. This method permitted reduc- 
ing the solution of the partial differential equations of the problem 
to the solution of linear differential equations of the second order 
with constant coefficients and reducing the determination of the con- 
stants of integration to the solution of an infinite completely regular 
system of linear equations. 

The obtained formulas determine the stiffness in torsion and bend- 
ing as a function of the geometric parameters of the section. 

At the same time, there are indicated the limits of applicability 
of the semiempirical formula of Bredt (reference 3) for the determina- 
tion of the stiffness in torsion of hollow thin- walled rods. 


1. TORSION OF A PRISMATIC ROD OF HOLLOW RECTANGULAR SECTION 

1. Statement of the problem . - The determination of the stress 
function U(x,y) for the torsion of a rod of a doubly connected cross 
section reduces, as is known, to the integration of Poisson's equation 


V 2 u = 


d 2 U 

c>x 2 



-2 


( 1 . 1 ) 


when the function U(x,y) becomes zero at the outer contour and assumes 
a constant value U 0 on the inner contour (fig. l). From the symmetry 

it is sufficient to find the function U(x,y) only for the part of the 
section ODEFBC . In order that the solution extend to the entire region 
of the cross section, it is necessary that on the lines DE and BC the 
normal derivatives of the function U(x,y) be equal to zero. 

_ ] izgib prismaticheskikh sterzhnei s polym pryamougol 1 nym 
•secheniem. Prikladnaya Matematika i Mekhanika, vol. 14, no. 3. 1950 
pp. 265-276. * * ’ 
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It is assumed that in the region OABC the function U(x,y) 
assumes the value U^(x,y) and in the region ODEG the value U 2 (x,y) . 

It is ideal to obtain the function U^(x,y) in the form 

Ui(x,y) = %(x,y) + ^(x.y) (i = 1,2) (1.2) 

where the functions ^(x^y) (i = 1,2) exists only in the region OAFG, 
l^foy) exists in the region OABC, and exists in the region ODEG. 

For the auxiliary functions ^ (x,y) and 4> i (x,y), the following 
equations are obtained: 

V 2 ^ = -2 V 2 ^ = 0 (i-1,2) (1.3) 

The following conditions must be satisfied: 

/a* V\ 

*i( x >o) = v^ry _ = + = 0 h( x ’ d i ) - u o d- 4 ) 

*2(°>y) = = *2 ( x j 0) + ^2( x >0) = ® ^2^2^ y ) = % ( 1,s ) 

The boundary conditions' for the determination of the functions 
^l(x,y) and g ( x , y ) are nonhomogeneous; however, following 

G. A. Grinberg (reference 4), set 


(x,y) = 



f k (x) sin * 2 ( x >y) = 2 v kd) sin 

k = 1 k= 1 


oo 

E 


kitx 


d c 


(1.6) 


For <I>- L (x,y) and $ 2 (x,y) the following conditions are then 
obtained: 




*i( d 2 ,y) = *2.(x,0) = l jr-±\ = 0 


dx 


x =d 2 


(1.7) 



M x > d i) = ) v k( d i) sin ITT “ u o 


k= 1 



MCA TM 1319 


3 


< M x > d i) = # 2(°>y) 



*2( d 2>y) 


> 

OP 

^ f k (d 2 ) sin _u 0 
k= 1 


( 1 . 8 ) 


It is ideal to obtain the functions $ 1 (x,y) and $ 2 (x,y) in the 

form 


oo 

^l(x,y) = <P k (x) sin ^> 2 (^y) = 



/ \ . kjrx 

w k (y) -Bin -=- 


k = 1 


k= 1 


(1.9) 


. Equations (l.3) to (l.9) completely determine the function U(x,y) 
in the region ODEFBC. 

2. Solution of the equations of the problem . - Making use of equa- 
tions (1.3) to (.1.9), the following equations are obtained: 


« / \ . , k?tx , kjtx / _ sk 2Uo 

f k (x) - A k sh dT - + B k ch dT ( _1 ) + 

1 1 


k:t 


t - <-> k ] 


4dn 


(kn) ; 


( 2 . 1 ) 


v k (y) = % ^ ♦ N k ch SI - <-l) k § ♦ [l - (-l) k ] ^ 


2 

(krt)8 


(2.2) 


» kW -%sh|S t c k c h ^ t (-l) k ^(l-=h^) 


(~1) ] 


OO 

•E 

P=1 


v p( d l) 


(pd-j/dg) 2 + k 2 


X 1 , knx , pjtx 


d 2 


p sh 


k sin 


d 2 


(2.3) 
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„ k(y) . I* sh & ♦ Fk oh ♦ (-!>* B! (l - ch ^ ♦ 


(-l) k | 




fp(i 2 ) 


J (pdg/d-L) 2 + k 2 

P= 1 


l 2 . kjry , . pjty 

r* p sh -~ L - k sin 


Then fp(dg) and Vp(d^) have the values 


p*a 2 ; _ pitd 2 , 2Ug + 4d 

” p + (p«) 


, v pjiao pita? 

f p (d 2 ) = Ap 8h -J- + Bp oh -J-S - (. 


, , PJtd-, pnd-i 2U n 4do^ r 

v p (d x ) , Mp sh + Up ch - 3 -i - (-i)P 5T + -ZV5 L 1 - ( ' 1)P J 


L 2 P * ( Plt ) 3 


(2.4) 


& [ p ' ^ 


(2.5) 


3. Determination of the- constants of integration . - For the func- 
tions f k (x), v k (y), <p k (x), and v k (y), the boundary conditions from 

equations (l.4), (l.5), (l.7), and (l,8) have the form 


f k '(a) = f k (0) + <p k (0) = 0 <P k (d 2 ) = <p k *(d 2 ) = 0 

(k = 1,2, . . . ) 

v k '(h) = v k (0) + v k (0) = 0 w k (d 1 ) = w k * (d x ) = 0 

(3.1) 

With the aid of these conditions the equations for the determina- 
tion of the constants of integration are obtained from equations (2.1) 
to (2.4). 


With the elimination of A k , C k , L^, and F k ' from the 

obtained equations, a set of two infinite systems of linear equations 
(reference 5) is obtained. The following notations are introduced: 


(-l)l 

% = S k d!d 2 - L ^ L 


knd«- 


sh 


N k = R k d 1 d 2 




kitdi 


sh 


(3.2) 


This set of two infinite systems of equations is then reduced to 
the form 
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oo 

= V V^P + 


^k — ^ ^p a kp + ^k B k ~ ^ ^ s p c kp + 7" k ( k 1*2, ... ) 

p=l (3.3) 


p=l 


where 


a kp * dl 2 p 2 + d2 2 k 2 


2k d l d 2 


. P nd l pnh P«(b" d i) 
sh — — sch ^ — ch = 


i f 2U 0 
= k [fV 


2 d 2 2 * 3 


+ 4 1 - (-l) k d l 


k d-i krtd 

csch 


kjld 2 2_ ( 2d l kltd 2^ 

d l ~ * 2 V krtd 2 2d l / 


(3.4) 


2k ' d l d 2 


, P* d 2 , pna , P«(a-d 2 ) 

c Iro = — — 5-5 5-5 sh — = — sch f — ch — -5 

* d 2 2 p 2 + d 1 2 k 2 d l d l d l 


_ lfl U 0 . 4 1- (-l) k d 2 , krtd l 2 L 2d 2 kjt M 

T * = k L^d! 2 * 3 k d l 08 d 2 ' ,2 V k ” d l w)_ 


(3.5) 


Systems (3.3) may he written in the form 


Zy ~ ^ ^ ^vpZp + (v — 1,2, ... ) (3.6) 

P =1 

where it is necessary to set 


Z 2n-1 = 

s k 

z 2n ~ R k 

A 2n,2m 

= 0 

A 2n,2m-1 “ 

c kp 

B 2n-1 = 

3k 

B 2n - r k 

A 2n-l,2m-l 

= 0 

A 2n-l,2m = 

a kp 


A few cases will he considered. 

1. The infinite system (3.6) is completely regular for b^d 2 and 
a^di since from equations (3.4), (3.5), and (3.7) 
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. k d 2 
C kp " it dp 


(pd 2 /di) 2 + k 2 


p = 1 


p = 1 


= ~ l cth 


d 2 kjtd-L 


l 2n-l,p 


a kp " Jt d 2 ) ' (p d:L /d 2 ) 2 + k 2 


p= 1 


P = 1 


If kltd 2 d l 
2 l cth d-L " kjtd 2 


h 


where the inequalities were used 

P* d l , pnb , P*( b " d l) _ , 

sh — ; — sch * — ch r * sh 


, P* d 2 , pita' , P*( a “ d 2) 
sh — : — sch *; — ch ; * sh 


PJtdl / 

f pnd-jj 

i 1 

(3.9) 

*3 “M 

, d 2 > 

n 2 

p«d 2 / 

' pjtd 2 ^ 


(3.10) 

A { 


h 

(o^ xs£ 00 ) 



(3.11) 


Hence, for any v 


a »p s i 


(3.12) 


2. The system (3.6) is regular for the particular cases where 
a>d 2 , b = dp and a = d 2 , b>dp (where a rectangular hollow section 

is present). 
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In this case 


]C Aan ' p = Cis< ” k a ' 11 , < pV ^) 2 + * 2 4 1 (3 ' 13) 

P= 1 P= 1 p = 1 


£W2><? »•« 

P= 1 P=1 p=l 

where the inequalities (3.1l) were used 

P nd l pnd-L pnCb-d-L)- pnd 2 , pjt(a-d 2 ) 

-T~ sch ch — T- <1 sh sch I— ch — < 1 


sh 


d 2 — d 2 - d 


2 ' 


d l ~“ d l 


(3.15) 


According to equations (3.13) and, (3.14) for any v 


Avp <1 (3.16) 

P = 1 

3. When a = b and d 2 = dp = d, a hollow square section (fig. 2), 

■a single infinite system of linear equations which is entirely regular 
for a/d ^ p. > 1, where |i is a finite number, is obtained in place of 
the systems (3.3). This results in 



where 


'k 


Y F p b kp + °k 

P = 1 


(k= 1,2, . . .) 


(3.17) 


(3.18) 
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’ 2U 0 

4 1 + 

(-l) k + 1 2 

(l - — 

th k*U 

[M 2 

kn3 

sh k« ~ n 2 

^ kn 

th *;] 


\\ ^ 2k 1 + e“ 2n (^“ 1 ) \ N 1 

2 _> T — 1 

p == 1 


1 + e" 2 *^" 1 ) 


(3.20) 


P=1 


where the inequalities (3.1l) were used and 

pna ^ pit(a-d) ^ 1 + e" 211 ^" 1 ) 

sh pit sch -- ch £ — 4 — L < 

d d 2 


(3.21) 


As an example, the case where a/d p = 3/2 is considered. From 
(3.20) is obtained 


go 

E 

P = 1 


b kp < 1 “ 0 = 0*5216 


9 = 0.4784 (3.22) 

The free term cc^ of the system (3.17) satisfies the inequality 


U, 


0 


|a k f$ 0.20264 - 0.06198 


(3.23) 


The values of the unknown F k with an excess is denoted by F k + 
and the values with a defect by F k “. 

Using the theory of regular and completely regular systems (refer- 
ence 5) and applying limiting values yields the following estimates for 

*k : 


U, 


0 


u, 


Ft” = 0.27547 - 5- - 0.09629 £ F-, $ 0.27588 - 0.09641 = F n + 

d 2 d 2 x 


u 0 - U n + 

F ' = 0.18234 - 0.11107 £ F ? $ 0.18300 - 0.11128 = F~ + 

2 d 2 2 d 2 2 


U, 


0 


V 


(3.24) 


F 3 " = 0.14762 -5- - 0.09749$ F 3 $ 0.14845 - 0.09775 = F* + 

0 d 2 ° d 2 5 


V U 0 + 

0.13440 - 0.09196$ F k “$ F k $ F k + $ 0.13739 -g- - 0.09288 


(k — 4,5, . . . ) 
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where Uq + is the value of Uq with an excess and Uq - is the value 
with a defect. 


4. Determination of the constant Uq. - 


is made of the theorem on the circulation of 
torsion (reference 6) 



For determining Uq, use 
the tangential stress in 

(4.1) 


where Cq is the inner contour FHMLF of the section (fig. l), 2q is 

the area enclosed hy this contour, G is the shear modulus, t is the 
angle of torsion per unit length, and T g is the projection of the 

tangential stress at any point of the contour Cq on the direction of 
the tangent to the contour. 

Substituting in (4.l) the value 

*.-(££-!£)<* < 4 - 2 > 


and maki ng a certain transformation, relation (4.1) is reduced to the 
form 



2(a-d 2 )(b-d 1 ) 


(4.3) 


Use of the obtained values of f^(x) and v-^(y) yields, from 
relations (1.2), (l.6), and (4.3), 


Uq = 


d l d 2 

dg (a-d 2 ) 4 d.^ (b-d-^ ) 


{ 


2ab - adn - bdo + 




kitd-^ 

sh —3 — sch 
d 2 


krtb 

d 2 


ch 


kn(b-d2_) 


+ 


k = 1 


(4.4) 



k = 1 


sch 


kna , kn(a-d 2 ) 


ch 


d. 



k = 1 


Pk + ^k I 
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where R k and S k are the constants of integration determined from 
the system (3.3); and 3 k and r k have the values (3.4) and (3.5). 
If the values of the coefficients R k and S k with excess and with 

defect are substituted in (4.4) and this equation solved for Uq, the 
values of Uq with excess and defect are obtained. 

5. Determination of the stress function . - According to (l.2) , 
(1.6), (1.9), and the obtained values of the functions f k (x), v k (y), 
<Pk(x), and w k (y) , the stress function has the form: 


U x (x,y) = 


■E 

P = 1 


v P (d i ) sin If csch sh + y( d i-y) + 


d l d 2 


V" ( _1 ) k IT sh sch ch sin 

/ , k d l d l d l d i 

k=l 


2U 0 'P (-l) k+1 . k *(x-d 2 ) . krty 

~ 2^ ch - d x sin + 

k =1 


for 0 ^ x ^ d- 


oo 

2 ii Y'' (-l) k sh kll(x " i 2 ) sin iss V s v p( d l> p(-l) p 

d 2 v — -i dl d l l—( (P d !/d 2 ) 2 + k 2 

iC — _L T) = 1 


Ui(x,y) = U 0 X + y(di-y) + 


(5.1) 

for dg ^ x ^ a 


d l d 2 


V (-l) k X sh -5^ sch If <* Kin I 22 


k = 1 


(5.2) 
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U 2 (x,y) = U 0 j- + x(d 2 -x) + 


\k *k 5^1 krtb ^ torfb-v) , k«x 


(-1) K X sh -gX sch ff fch 


U 2 (x,y) = 


fp( d 2 ) sin csch ~T^ Sh d^ + x ^ d 2“ x ) + 
U- 1 1 1 


P =1 


(-l) k ^ sh ^ sch £2* ch Sin *52 + 

k d 2 d 2 d 2 d 2 


L )k+1 , k "(y- d l) . knx 

-i- ch sin - — + 

k d 2 d 2 


2^2^ (_i)k ah kff (y- d l) gin knx^ f p ( d 2 )(-l) P P 


P-1 


(p d 2 / d l) 2 + k 2 


where f (d 2 ) and v p (d-j_) have the values 


(_i]p+i 


f P< d 2> - ^ 


2U n 4d n 2 r - Ll "l 

+ £ [* + . <- 1>P+1 ] - 


Sp d i d 2 sh -j- sch 


(5.5) 
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... . Prtd l „ P*b _ P«( b - d 

R p d l d 2 sh Sch d^T Ch d^ 




(5.6) 


By the substitution of the obtained values of the stress function 
U(x,y) in the general formulas, the stress and stiffness in torsion are 
obtained. 


As an example the torsional stiffness and stress of a rod with a 
hollow square section (fig. 2) will be determined. 

6. Determination of the torsional stiffness . - The substitution 
in equations (5.1) to (5.6) of b = a and d-j_ = d 2 = d and the use 

of the formula for the stiffness of a section 


d d 


C = 8G (a-d) 2 u 0 + 1 

j U(x,y) dxdy + 2 J 

d y j 

U(x,y) dx 

L J o 

U J 

0 ( 

5 c 

i (6. 1 ;. 


yields, after integration, 
C = 8Gd 4 l 


( a . a ^0 + 1 ( a . -A 

U ) d d 2 3^d ) 


( 6 . 2 ) 


7 £ ?[ 1 -^ th T s -s 8hM ^ lsoh ¥! + 


k- 1 , 3 , o 


4 


oo 

p [l + sh k„ soh sh idazil . 


k=l,3, . . 


seh ch Mali 

d d 




p(-l)P +1 f p( d ) _ 


k 2 d 2 


p =1 

(continued on the following page) 
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p — 1 ^ 3 j « • • 


where Uq is determined from the equation 



the unknown constants F^ being determined by inequalities (3.24) and 
a k having the value (3.19). 


Substitution of the obtained values of the coefficients F k + and 
F k " in (6.2) yields the upper and lower limits of stiffness. The 
coefficients with a defect F k ~ will correspond with the lower stiff- 
ness limit C~ and the coefficients with an excess F k + will corre- 
spond with the upper stiffness limit C + . 
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TABLE I 


a/d 

c 0 * 

co + 

Co" 

A 

5 

1.5 

8.0 

11.052 

11.051 

38.4 

0.009 

2.0 

27.0 

32.952 

32.949 

22.0 

.008 

2.5 

64.0 

73.780 

73.775 

.15.3 

.008 

3.0 

125.0 

139.518 

139.509 

11.6 

.007 

3.5 

216.0 

236.164 

236.150 

9.3 

.007 

4.0 

343.0 

369.716 

369.697 

7.8 

.005 

5.0 

729.0 

771.542 

771.509 

5.8 

.004 

10.0 

6859.0 

7035.066 

7034.902 

2.6 

.002 

20.0 

59319.0 

60034.051 

60033.323 

1.2 

.001 


In table 1 are given the relative values of the stiffness computed 
by formula (6.2): 

C Q + = C + /Gd 4 C 0 " = C“/Gd 4 
and the maximum relative error 

5 S! (c + - C + )/C" 

in percent for different ratios a/d. 

For comparison there are also given in table 1 the values of the 
relative stiffness computed by the semi-empirical formula of Bredt 
(reference 3). 



and the different A = (C + - C*)/c + in percent. 

From table 1 it is seen that the semi-empirical formula of Bredt 
(6.5) gives sufficiently close results only for thin-walled rods for 
which a/d 5. For thick- walled rods, however, for which a/d < 5 
the Bredt formula is not applicable; for a/d = 4 it gives an error 
of 8 percent which rapidly increases. 

7. Determination of the stresses . - With the use of expressions 
(5.l) to (5.2) for the stress function, the stresses are readily 
obtained by the usual formulas of the theory of elasticity. 

Leaving out the computations gives the following expressions for 
the stresses at the points (a,0) and (a,d) for the rod with a hol- 
low square section represented in figure 2: 
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X z (a,0) = 


k = 1 


l) k F k sch ~ sh kir^ Grd (7.l) 


■} 


X z (a,d) =J - 4 ; - 1 + it ^ ^ F k sch -^p sh knl 

*“ fcl J 


Y z (a,0) = Y z (a,d) = 0 


Ord 


(7.2) 


(7.3) 


After the substitution of the coefficients F k + and F k ", the 
upper and lower limits of the stresses X z (a,0) and X z (a,d), where 
the upper stress limit X z + will correspond with the coefficients with 
an excess F k + (k=l,2, . . .), are determined. 


In table 2 are given the computed values of the stresses X z (a,0) 
and X z (a,d). In the first column of each stress are given the values 
with an excess and in the second column, with a defect. 


TABLE 2 


a 

d 

X z (a,d) 

Grd 

X z (a,0) 

Gxd 

Xz*max 

Gt& 

1.5 

0.3149 

0.3143 

2.0212 

2.0212 

2.0 

2.0 

0.7325 

0.7322 

2.6326 

2.6324 

2.25 

2.5 

1.2085 

1.2082 

3.1798 

3.1797 

2.666/ 

3.0 

1.7050 

1.7048 

3.6974 

3.6972 

3.125 

3.5 

2.2068 

2.2066 

4.2048 

4.2046 

3.6 

4.0 

2.7092 

2.7089 

4.7088 

4.7085 

4.0833 

5.0 

3.7131 

3.7129 

5.7131 

5.7129 

5.0625 

10.0 

8.7204 

8.7201 

10.7204 

10.7201 

10.0278 

20.0 

18.7236 

18.7234 

20.7236 

20.7234 

20.0131 


In the same table the values of the maximum stresses X z *, 
computed by the formula of Bredt, are given for comparison. . 


k z max* 




(7.4) 
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From table 2 it is seen that in the determination of the maximum 
stresses, the approximate formula of Bredt may be applied only for very 
thin- walled hollow rods. 

When a/d = 20, the Bredt formula (7.4) gives an error of 
3.5 percent. 

With decrease in the ratio a/d this difference increases; for 
a/d = 10 it is equal to 7 percent, for a/d = 5 it is equal to 
13 percent, and for a/d = 3 it is equal to 18 percent, etc. 


II. BENDING OF A PRISMATIC ROD OF HOLLOW RECTANGULAR SECTION 

8. Statement of the problem . - The stress function F(x,y) in 
bending, as is known, satisfies the equation 


V 2 F = 


d 2 F 

dx 2 


d 2 F 

dy 2 


Pv 

l(l+v) 


(y-yo) - 2i f '(y) 


( 8 . 1 ) 


within the region of the section and the condition 

g = [x2 - 2xxo - f (y)] g (8.2) 

on the contour of the section where P is the bending force applied to 
the free end of the rod at the center of gravity of the section, v is 
the Poisson coefficient, Xq, Jq are the coordinates of the center of 

gravity of the section, I is the axial moment of inertia of the sec- 
tion about the y axis; the arbitrary function f(y) is to be deter- 
mined from the conditions at the contour. 


On account of the symmetry (fig. 3) it is sufficient to find the 
function F(x,y) only for the part ODEFBC of the section. 

In order to extend the solution over the entire region of the 
cross section, it is required, on the basis of the membrane analogy 
(reference 7), that the function F(x,y) becomes zero on the vertical 
axis of symmetry and the derivative c)F/dx becomes zero along the 
horizontal axis of symmetry. 

By condition (8.2) on the section contour, the function F(x,y) 
is determined with an accuracy up to a constant term. For the cross 
section which is a doubly connected region, the number of constant 
terms is equal to two; and for their determination use is made of the 
theorem on the circulation of the tangential stresses in bending 
(reference 6). 
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It is assumed that in the region OABC the function F(x,y) 
assumes the value F d (x,y) and in the region ODEG the value F 2 (x,y). 

It is ideal to obtain the function F d (x,y) (i = 1,2) in the form 

*i(x,y) = ¥ ± (x,y) + <t>i(x,y) (i = l,2) (8.3) 

where the functions %(x,y) (i = 1,2) exist only in the region OAFG, 

the functions Y x (x,y) exist in the region OABC, and ¥?(x,y) • in the 
region ODEG. ^ 


For the auxiliary functions ¥ d (x,y) and ^tay) (i = 1,2), 
setting f(y) = 0, the following equations are obtained: 

V 2 Y i= K(y-b) V 2 $i = 0 (i = 1,2) ( K = f (8 ‘ 4 > 

where the following conditions must be satisfied: 


»l(x,0) 



^i(o,y) 


+ $l(o,y) 


= 0 


Yl(x,di) = C 1 


(8.5) 


l 2 (x,b) = Y g (0,y) = ^ 2 (x,0) + ^2 ( x ^°) = 0 
y 2( d 2>y) = Jv (2 a -d 2 )(b-y)d 2 


(8.6) 


By the method described in the first section set 


°° CO 

f i( x >y) = sin ^ M x >y) = 

k=l k=l 


v k(y) 


knx 

SID 


(8.7) 


Then 


♦lfey) = v^r/ x=d2 = *i( x >°) = M x > d i) = ) v k (d-L) 


knx 

d 2 


sin - C x 


k = 1 




: (x ; d l) = - *z(o,y) - 


( 8 . 8 ) 


0 
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f k( d 2) sin ■ 4 i ( 2a - d 2)( b -y) d 2 (8.9) 


z_^ 

k = 1 


form 


It is ideal to obtain the functions ^(x,y) and ^tay) in the 


*l(x,y) = $ k (x) sin 

k = 1 


^>2C x ^y) - ^ v k (y) sin 

k=l (8.10) 


Equations (8.4) to (8.10) determine F(x,y) in the region ODEFBC. 


The constant C-^ is determined with the aid of the theorem on the 

circulation of the tangential stresses in bending for the inner contour 
of the section: 



9. Determination of the bending stresses . - The solution of the 
equations of this problem is analogous to the solution of the equations 
of the problem of the torsion of a rod. 


With the omission of the computations the values of the obtained 
stresses axe: 

for the region OAFG 


X z (x,y) = ^ (2a-x)x + ^ (Sy 2 ^ 2 ) - £ (Zy-d-^b + 


2^ P v 


09 

Jt \ . _ . kita . km(a-x) kmy 

‘dj ) “k SCh df Ch df* C0S df 


L^r 

k = 1 


00 

Jt V A , ( a \ ' kxx kjty 

di k( l} Sin if CSCh “df ch df - 


k = 1 


(continued on following page) 
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k = 1 

2. V' k(-l) k sh kir(x - d ^ cos *2 V T P (dl) P( ' 1)P 

d 2 d l d l l_j (pd-L/dg) 2 + k 2 

k= 1 P= 1 


sch |= sh ^ 8in ^ . 

k = 1 


OO 

p= 1 


pv_(dT) cos csh ^---- sh £22 - 
y dg dg d2 


2. V" t , 1)k oh k »<*- d z> sin is ^ V d i> + 

d 2 A d l d l L (pll / d 2 ) 2 + k 2 

p=l k= 1 


3 <-l>* sh k ’ ,(X - d2 ) - sis 


2Ci 

id L 

k= 1 


sin 


for the region GFBC 


X z( x ^) “ 21 ( 2a_x ) x + + 61 1 + v " d l ^ 

TTT (2y-dx)b + i Y “k soh If oh cos ff 


(9.1) 


_P_ 

21 


k = 1 
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Y z (x,y) = 




i2S sh sin |S (9.2) 


k = 1 


for the region APED 


X z (x,y) = -Jj (2a-x)x - 

k =1 


k d 2 b 


kBk + (-D* -f- 


I (2a-d 2 ) + 


2P v 


I 1 + v 


(kn) 


2 (1 + (-l) k+1 ) 


knb , kn(b-y) . knx 
- — - ch — sin _ — _ 


> csch -z — - ch 


it ( 2a " d 2) x - TTT (d 2- x ) x 


Y z (x,y) = - jkDjj. + (-l) k ^ 


2P v 


d 2 


I 1 + v 




k= 1 

_ 

(i + (-i) k+ i) 


| (2a-d 2 ) + 


! , knb , kn(b-y) knx 

> csch -t— sh — ---■ ■■ < cos — - 


(2a-d 2 ) (b-y) - ^TT ( b “ y) ( d 2 _2x ) 


L 2 

(9.3) 


In these relations the constants of integration and are 

determined from the infinite completely regular system of linear 
equations. 


Translated by S. Reiss 
National Advisory Committee 
for Aeronautics 
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